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We use the exact strong-interaction limit of the Hohenberg-Kohn energy density functional to 
approximate the exchange-correlation energy of the restricted Kohn-Sham scheme. Our approxima- 
tion corresponds to a highly non-local density functional whose functional derivative can be easily 
constructed, thus transforming exactly, in a physically transparent way, an important part of the 
electron-electron interaction into an effective local one-body potential. We test our approach on 
quasi-one-dimensional systems, showing that it captures essential features of strong correlation that 
restricted Kohn-Sham calculations using the currently available approximations cannot describe. 



In principle, Kohn-Sham (KS) density functional the- 
ory (DFT) [I] [5] should yield the exact ground-state den- 
sity and energy of any many-electron system, including 
physical situations in which electronic correlation is very 
strong, representing them in terms of non-interacting 
electrons. Currently available approximations for KS 
DFT, however, fail at properly describing systems ap- 
proaching the Mott insulating regime, the breaking of 
the chemical bond, and localization in low-density nan- 
odevices, to name a few examples (3] 0]. Artificially 
breaking the spin (or other) symmetry can mimic some 
(but not all) strong-correlation effects, at the price of 
a wrong characterization of several properties and of a 
partial loosening of the rigorous KS DFT framework. 

Indeed, it is very counterintuitive that strongly- 
correlated systems, in which the electron-electron repul- 
sion plays a prominent role, can be exactly represented 
in terms of non-interacting electrons. For this reason, 
several authors [5H2"0] have used accurate many-body 
solutions of prototypical strongly-correlated systems to 
obtain (by inversion) and characterize the exact non- 
interacting KS system. The exact properties needed to 
describe strong correlation in KS DFT have also been set 
in a transparent framework [3TJ [55] . These works made 
it all the more evident how difficult it is to find adequate 
approximations of the exact KS system [4 1, so that, albeit 
theoretically possible, it may seem unrealistic to describe 
strongly-correlated systems with KS DFT. 

Here, we address this skepticism by showing that the 
strong-interaction limit of the Hohenberg-Kohn (HK) en- 
ergy density functional yields approximations capturing 
strong-correlation effects within the non-interacting re- 
stricted self-consistent KS scheme. 

The Letter is organized as follows. First, we intro- 
duce the formalism, using the strong-interaction limit 
of the HK functional to transform exactly an important 
part of the many-body interaction into an effective lo- 
cal one-body potential, in a physically transparent way. 
We then present pilot self-consistent Kohn-Sham calcu- 
lations, showing that this potential is indeed able to cap- 
ture strong-correlation effects way beyond the reach of 
present KS DFT approximations. As a prototypical ex- 



ample, we look at the "2kp — > 4fcp" crossover of electrons 
confined in quasi-one dimension (Q1D). This crossover is 
entirely due to the dominant particle-particle repulsion 
that tends to localize the charge density, destroying the 
non-interacting shell structure, as it happens in many 
strong-correlation phenomena. The interest of these re- 
sults goes beyond quasi-one-dimensional systems, be- 
cause the latter are a valid test lab for three-dimensional 
DFT, as clearly discussed in [23]. Our approximation 
turns out to be qualitatively right, and quantitatively 
very accurate for the ionization energies, although less ac- 
curate for the ground-state density. We thus conclude by 
discussing the inclusion of higher-order corrections and 
strategies for extending the self-consistent calculations 
to two and three dimensions. Hartree (effective) atomic 
units are used throughout. 

Strong-interaction limit- In Hohenberg and Kohn's for- 
mulation [T] the ground-state density and energy of a 
many-electron system are obtained by minimizing with 
respect to the density p(r) the energy density functional 

E[p] =F[p}+ J drv cxt (r)p(r), (1) 

where w ox t(r) is the external potential and F[p] is a uni- 
versal functional of the density, defined as the minimum 
of the internal energy (kinetic energy T plus electron- 
electron repulsion V ee ) with respect to all the fermionic 
wave functions "J that yield the density p(r) [54] , 

F[p] = min(*|f' + V r ee |*). (2) 

To capture the fermionic nature of the electronic density, 
Kohn and Sham [5J introduced the functional T s [p] by 
minimizing the expectation value of T alone over all the 
fermionic wave functions yielding the given p [24], 

T s [p] = min(tf |f|*>, (3) 

thus introducing a reference non-interacting system with 
the same density of the physical, interacting, one. The re- 
maining parts of F[p], defining the Hartree and exchange- 
correlation functional, i?Hxc[p] = F[p]—T s [p], are approx- 
imated. 
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The strong-interaction limit of F[p] is given by the 
functional V e s e CE [p], defined as [25H2E] 



min(*|y ee |*), 



(4) 



where the acronym "SCE" stands for "strictly-correlated 
electrons" [25]. The functional V e S e CE [p] is the minimum 
of the electronic interaction alone over all the wave func- 
tions yielding the given density. It has been first studied 
in the seminal work of Seidl and coworkers [25U27] . and 
later formalized and evaluated exactly in Refs. [28H3T1 

More recently, it has been suggested that a "SCE 
DFT", in which the functional F[p] is decomposed as 
F[p] = Vl CE [p] + E kd [p] HE [31 133], and t fl e so- 
called kinetic-decorrelation energy Ekd[p] is approxi- 
mated, could be a complementary alternative to KS 
DFT for systems in which the electron-electron repulsion 
largely dominates over the electronic kinetic energy. In- 
deed, SCE DFT works well for low-density many-particle 
scenarios [291 132] , but it requires that one knows a priori 
that the system is in the strong-interaction regime, and it 
fails when the fermionic shell structure plays a role [2*5] . 
It also misses several appealing features of KS DFT, e.g. 
the possibility to yield (at least in principle) the exact 
ionization energy. More generally, orbitals and orbital 
energies, crucial for chemistry and solid state physics, 
are totally absent in SCE DFT (while captured by KS 
DFT). 

SCE as a functional for KS DFT - To combine the 
advantages of KS DFT and SCE DFT, here we use the 
functional V e g CE [p] to approximate -Eh~xc[p], 



#HxcM 



V. 



SCE 



[p\- 



(5) 



Equation ([5| amounts to approximating the constrained 
minimization over '5 in the HK functional ^ with the 
sum of two constrained minima, 



mm(V\T + V ee \V) 

<rp 



min(*|T|*) 



-min(*|14e|*). (6) 

>p 



This new "KS SCE" approach treats both the kinetic en- 
ergy and the electron-electron repulsion on the same foot- 
ing. Standard KS DFT emphasizes the non-interacting 
shell structure, treated accurately through the func- 
tional T s [p], but it misses the features of strong corre- 
lation. SCE DFT is biased towards localized "Wigner- 
like" structures in the density, missing the fermionic shell 
structure. Many interesting systems lie in between these 
two limits, and their complex behavior arises precisely 
from the competition between the fermionic structure 
embodied in the kinetic energy and correlation effects 
due to the electron-electron repulsion. By letting these 
factors compete in a self-consistent KS procedure, one 
might be able to get at least a qualitative description of 
several complex phenomena, amenable to improvement 
by corrections in the same spirit of standard KS DFT. 



General features of KS SCE - First, notice that for a 
given density p, the right-hand side of Eq. |6]) is always 
less or equal than the left-hand side. Even if minimizing 
our energy functional with respect to the density will not 
yield the exact p [as Eq. ^ is an approximation], it is 
easy to prove that our final total energy is a lower bound 
to the exact one. 

From the scaling properties [34] of F[p], T s [p] and 
Ke CE [/°] h derives that the approximation of Eq. ^ is 
accurate both in the weak- and in the strong-interaction 
limits, while probably less precise in between. By defin- 
ing, for electrons in D dimensions, p 7 (r) = r y D p(jr), 
where 7 > 0, we have T s [p 7 ] = 7 2 T s [p], F e s e CE [p 7 ] = 
7 K S e CE [p] M, and F[p 7 ] = 7 2 ^i/>], where F lh [p] 
means that the Coulomb coupling constant in F[p] is 
rescaled with I/7. It is then easy to see that both sides of 
Eq. ([6| tend to T s [p 7 ] when 7 — >■ 00 (high density or weak 
interaction) and to V e | CE [p 7 ] when 7 — » (low density or 
strong interaction). 

Since KS SCE tends to the exact density and energy in 
the strong-interaction limit, the corresponding KS poten- 
tial should have the features that are expected for a KS 
description of strong correlation (6[ H3J [16] . We discuss 
first why, physically and mathematically, the SCE po- 
tential [Eqs. ([7])- ^ below] is expected to have these fea- 
tures, which we then test practically with self-consistent 
calculations in Q1D. 

Physically, the functional [p] portrays the strict 

correlation regime, where the position r of one electron 
determines all the other N — 1 electronic positions 17 
through the so-called co-motion functions, 17 = f; [p] (r) , 
some non-local functionals of the density [351 [3T] [32 [35] . 
Therefore, the net repulsion on an electron at position r 
due to the other N — 1 electrons depends on r alone. Its 
effect can then be exactly represented [551 1311 132] by a 
local one-body potential [15] . 



V'CscE[ ( o](r) 



JV 

E 

i=2 



fJp](r)P 



(7) 



The physical meaning of Eq. is very transparent: at 
each position r, V«sce [p] ( r ) exerts the same force as the 
net electron-electron repulsion. We also have [3_U [31] 



SV^[P] 
Sp(r) 



«SCE[p](r), 



(8) 



so that the approximation of Eq. ^ corresponds to 
model the exchange-correlation potential v xc [p] (r) of KS 
DFT as v xc [p](r) sa u SCE [p] (r) - v H [p] (r) , where v H [p](r) 
is the Hartree potential. The functional V e g CE [p], being 
essentially a classical repulsion energy, favors localized 
charge densities. When evaluated with a delocalized den- 
sity p, its functional derivative ^ as a function of r dis- 
plays strong variations pushing electrons towards local- 
ization. Otherwise said, Eqs. ^ transfer the effects 
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FIG. 1: (color online) Self-consistent densities for N — 2 elec- 
trons in Q1D [hamiltonian of Eqs. (|9])-(10l with b = 0.1 and 



«oxt(a:) = sw i ], in units of the effective confinement length 
L — 2lu~ 1 / 2 (here and in the following figures). The exact 
results are compared with KS LDA and KS SCE approxima- 
tions. At large L the KS LDA calculations do not converge, 
while KS SCE approaches the exact solution. 



of strong-correlation into a physically meaningful, effec- 
tive local potential, expressed as the functional derivative 
of a rigorous KS density functional. 

While KS SCE does not use explicitly the Hartree func- 
tional, the correct electrostatics is still captured, since 
V e g CE [ j o] is the classical electrostatic minimum in the 
given density p. Moreover, the potential «scE[p](r) stems 
from a wave function (the SCE one [23 [35]) and is there- 
fore completely self-interaction free. Similarly, we expect 
-CgcE [p] (r) to have a derivative discontinuity that will be 
analyzed elsewhere [55]. 

Self- consistent KS SCE calculations in Q1D - As a 
pilot test of the approximation of Eq. ^ , we consider N 
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FIG. 2: (color online) Self-consistent exchange-correlation po- 
tentials for the same system of Fig. [I] For clarity, the poten- 
tials for L — 29 have been shifted by —4. 



electrons in a thin quantum wire with hamiltonian 
1 N 



H 



dxf 



N-l N N 

^2 X! w b(\ x t~ x j\) + ^2 v cxt(xi), 
i—1 1 i—1 i — 1 

. (9) 

where the effective Q1D interaction is obtained by inte- 
grating the Coulomb repulsion on the lateral degrees of 
freedom [37], 



(10) 



The parameter b fixes the thickness of the wire, and 
erfc(x) is the complementary error function. The inter- 
action wt>{x) has a coulombic tail, Wb(x — > oo) = 1/x, 
and is finite at the origin, where it has a cusp. 

The co-motion functions fi(x) for N electrons can be 
constructed from the density p(x) [25], [30l [31] : 



M*) = 



N- 1 [N e (x)+i-l] 



x > ajv+i-i, 

(11) 



N- 1 [N e (x)+i- 1 - N 
where the function N e (x) is 

N e {x) = [ p{x')dx', (12) 

J — oo 

and afc = N~ 1 (k). Equation Q becomes in this case 

JV 

v'sCe[p\( x ) = ^2w' b (\x - fi(x)\)sgn{x - fi(x)). ( 13 ) 

i=2 

We then solve self-consistently the restricted KS equa- 
tions in the KS potential vks(%) = v ex t(x) + vsce[p]{x), 
where vsce[p\(x) is obtained by integrating Eq. (13 1 with 
the boundary condition wsce[p](M °°) = 0. 

Here, we aim at showing that this KS SCE approach 
captures essential features of strong correlation out of 
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FIG. 3: (color online) Self-consistent KS SCE densities for 
N = 4 electrons in Q1D [hamiltonian of Eqs. |9]|-(10l with 
b = 0.1 and v cxt {x) = \uj 2 x 2 \. 



FIG. 4: (color online) The KS eigenvalue obtained in the self- 
consistent KS LDA and KS SCE calculations for the same 
N — 2 electron problem considered in Fig. [TJ plotted against 
the negative of the exact ionization energy, En — En-i- 



reach for standard restricted KS calculations. A simple 
but very representative example is provided by Abed- 
inpour et. al. |38j . who considered the external har- 
monic confinement v ex t(x) — ^lu 2 x 2 , and performed self- 
consistent KS calculations within the local density ap- 
proximation (LDA, [39]). Fig. [TJ shows our results for 
N = 2, together with accurate exact values [38]: as ex- 
pected, KS LDA works well when correlation is weak or 
moderate, a case characterized by relatively small val- 
ues of the effective confinement length L = 2lo~ 1 / 2 . As 
correlation becomes stronger (large L), KS LDA cannot 
describe the "2fcp — > Akp" crossover, simply reflected 
by the doubling of the number of peaks in the density. 
Indeed, a local or semilocal functional of the density can- 
not describe this crossover [38] , and exact exchange per- 
forms even worse. To localize the charge density, the self- 
consistent KS potential must build a "bump" (or barrier) 
between the electrons [35]. This "bump" was discussed 
in Refs. [TH] and [5] it is expected to be the key feature 
enabling a KS DFT description of the Mott transition 
and the breaking of the chemical bond, and it must be a 
very non-local effect [6] . We see in Fig. [TJ that the self- 
consistent KS SCE densities, although, as expected, less 
accurate in between the weak- and the strong-interaction 
cases, capture the transition to the strongly-correlated 
regime, thus building, at least partially, the "bump" in 
the self-consistent KS potential. This is confirmed by the 
exchange-correlation potentials reported in Fig. [2] we see 
that the "bump" is clearly present and gets closer to the 
exact one as the strong-interaction regime is approached. 
The long-range part of the SCE potential is also remark- 
ably accurate, as expected from the fact that the SCE 
functional is self-interaction free. Fig. [3] displays the KS 
SCE densities for N = 4 electrons: again, we clearly see 
the crossover from two peaks (the non-interacting shell 
structure) to four peaks (charge localization). 

In the exact KS theory, the highest occupied KS eigen- 
value is equal to minus the exact ionization potential 



I = E N -! - E N 001 ETJ. In Fig. g we plot the KS 
LDA and KS SCE eigenvalues for N = 2, as a function 
of the exact difference — Eff—\ for several harmonic 
confinement strengths. We see that KS SCE is remark- 
ably accurate [4*4"] . 

Conclusions and Perspectives - The exact strong- 
interaction limit has the promise of extending KS DFT 
applicability to strongly-correlated systems, while retain- 
ing the appealing properties of the Kohn-Sham approach. 
In Q1D the computational cost of KS SCE compares to 
KS LDA. Crucial for future applications is calculating 
Ke CE [/°] an d £sce[p](i*) also for general two- and three- 
dimensional systems. An enticing route towards this goal 
involves the mass-transportation-theory reformulation of 
the SCE functional [3Tj, in which V e s e CE [p] is given by the 
maximum of the Kantorovich dual problem, 



u(r)p(r)dr 



N JV-1 N ) 

i=l t=l j>i 1 1 Jl ) 



where u(r) = CgcE[/°]( r ) + C, and C is a constant [31] . 
This is a maximization under linear constraints that 
yields in one shot the functional and its functional deriva- 
tive, and can also inspire approximate and simplified ap- 
proaches to the construction of Vj|, CE [/9] and vscE[p]( r ), a 
critical step for the computational cost of KS SCE. An- 
other important issue is to add corrections to Eq. ([5]). 
One can, more generally, decompose F[p] as 



F[p] = T s [p] + V^[p] + T c [p] + V e d e [p] 



(14) 



where T c [p] (kinetic correlation energy) is the difference 
between the true kinetic energy and T s [p], and V^ e [p] 
(electron-electron decorrelation energy) is the difference 
between the true expectation of V ee and V^I, CE [p]. A 
"first-order" approximation for T c [p] + V^ e [p] can be, in 
principle, included exactly using the formalism developed 
in Ref. 1351 but other approximations, e.g. in the spirit of 
Ref. H21 can also be constructed. 
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